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Spurious 'active longitudes' in parametric models of heavily spotted 
eclipsing binaries 



i 



in 
o 
o 

(N 

u 
Oh! 

< 
oo 



> 

m 

o 
in 
o 

Of 

6 



S.V. Jeffers 

School of Physics and Astronomy, University of St Andrews, North Haugh, St Andrews, Fife KY16 9SS, U.K. 



ABSTRACT 

In this paper, size distributions of starspots extrapolated from the case of the Sun, are mod- 
elled on the eclipsing binary SV Cam to synthesise images of stellar photospheres with high 
spot filling factors. These spot distributions pepper the primary's surface with spots, many 
of which are below the resolution capabilities of eclipse mapping and Doppler imaging tech- 
niques. The lightcurves resulting from these modelled distributions are used to determine the 
limitations of image reconstruction from photometric data. Surface brightness distributions 
reconstructed from these lightcurves show distinctive spots on the primary star at its quadra- 
ture points. It is concluded that two-spot modelling or chi-squared minimisation techniques 
are more susceptible to spurious structures being generated by systematic errors, arising from 
incorrect assumptions about photospheric surface brightness, than simple Fourier analysis of 
the light-curves. 
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1 INTRODUCTION 

Doppler imaging has shown that rapidly rotating RS CVn binary 
systems, such as SV Cam (F9V + K4V, P rot =0.59d) frequently 
show spots at high and polar latitudes. A successful theoretical 
model of the formation of polar spots assumes that magnetic flux 
from decayi ng active regions is sw ept toward the poles by merid- 
ional flows (Schriiver & Title 2001). However, in order to produce 
polar spots bipolar active regions have to emerge at a rate approx- 
imately 30 times faster than in the case of the Sun, implying that 
the photospheres of active stars should be peppered with a large 
number of small spots. 

Jeffers et al.(2005) used spectrophotometric data from the 
Space Telescope Imaging Spectrograph on board the Hubble Space 
Telescope to eclipse map the inner face of the primary of the RS 
CVn SV Cam. These observations and the HIPPARCOS parallax 
showed that the surface flux in the eclipsed low latitude region 
is approximately 30% lower than computed from the best fitting 
PHOENIX Allardet al. (2000) model atmosphere. This flux deficit 
can only be accounted for if approximately 28% of the primary's 
surface is peppered with dark spots too small to be resolved through 
eclipse-mapping techniques where the resolution limit is approxi- 
mately 1.5 degrees. 

This paper extends the work of Jeffers et al. (2005) by ap- 
plying the extrapolated solar spot size distribution of Solanki 
1 1999) to a hypothetically immaculate primary star of SV Cam. 
We determine how different degrees of spot coverage influence the 
shape of the binary system's lightcurve, and how accurately these 
lightcurves are reconstructed into surface brightness distributions 



using the Maximum Entropy (Max Ent) eclipse mapping technique 
I Collie r Camera J 199llCoriier Cameron & Hilditohl 19971) . 



2 MODELLED DATA 

2.1 Size distribution of spots on active stars 

The variable nature of the spot coverage of active stars makes the 
quantification of their spot size distribution an intriguing problem. 
In the ca se of the Sun, the sp ot size distribution has been deter- 
mined by Bogda n~t alJll988l) from direct observations taken from 
the Mount Wilson white-light plate collection covering the period 
1917-1982. When plotted on a log-log scale the size distribution 
of sunspots is parabolic, implying a two-parameter log-normal dis- 
tribution, which cannot be represented b y a single parameter dis- 
tribution such as a power law. Following Bogdar Tet alJ fl988) the 
number of sunspots, N, as a function of the solar surface area, A, is 
equated as, 



dN _ (dN\ 



exp 



(}nA-ln{A))' i 
2 In a a 



(1) 



where the constants (A) and a a are the mean and geometric stan- 
dard deviation of the log-normal distribution, and ( 4¥) the 
maximum value reached by the distribution. For the case of the 
Sun, these values are tabulated in TableQwhere set 1 is for an in- 
active Sun and set 2 is for an active Sun. 

The extrapolation of the solar spot size distribution to active 
stars can be represented by different extrapolations e.g. [golanki 



2 S.V.Jeffers 



Set 


1 


2 


3 


4 


5 


6 


7 


Spot Coverage 




0.3% 


1.6% 


6.1% 


18% 


48% 


100% 


& A 


3.8 


5.0 


6.8 


9.2 


12.2 


15.8 


20.0 


(dN/dA) miM 


5 


25 


65 


125 


205 


305 


425 


Aspot.tot (XlO -6 ) 


320 


3000 


16000 


61000 


186000 


484000 


10 6 




200 


400 


1000 


3000 


6000 


10000 


10000 


Solar: 


min 


max 













Table 1. Tabulation of the input parameters to th e log-normal siz e distribution of star spots, equation The parameters derived bv lBogdanet'aii n'988) for 
the Sun are data sets 1&2, and those calculated bv lSolank] 1 19991) for active stars are data sets 3-7. 



( 1999) and Sola nki & Unruhl <2004 . As long as the size distribu- 
tion shows many sub-resolution spot that are not clumped together, 
the exact representation of the extrapolat ed size distribu tion is of 
minor importance. The extrapolation of Solankj) i 199911 is used, 
which is determined by analysing how the observable parameters 
of the log normal distribution, (dN/dA) maa; , (A), and a a, change 
with increasing magnetic activity. The extrapolated size distribu- 
tions cover a range of degrees of spottedness, from the one at solar 
mini mum to a hypot hetically completely spotted star. As adopted 
in the lSolankil (1999) model, the penumbral diameter is included in 
the spot area (in the umbral to penumbral area ratio of 1:3) and is 
assumed to be independent of spot size. In all cases the minimum 
spot area is 1.5x 1O _6 A0, where Aq is the surface area of the visi- 
ble solar hemisp here. The input param eters for the lo g-normal dis- 
tributions of Bo gdanetalJ ( 1988) and lSolankil fl999) are tabulated 
in Table Q sets 1 and 2 are for the inactive and active Sun, and 
sets 3 to 7 correspond to active stars with a spot coverage of 1.8%, 
6.1%, 18%, 48% and 100% of their surface. 



2.2 Generation of randomly distributed spots 

The binary eclipse mapping code DoTS iCollier Cameronlll997t) 
was used to synthesise spot maps which follow a log-normal size 
distributions on the surface of an immaculate SV Cam's primary 
(primary: Ri=1.24R , T e// =6038K, secondary: R 2 =0.79R Q , 
T e //=4804 K) to incorporate the surface geometry and radial ve- 
locity variations of tidally distorted close binary stars. The input 
parameters to DoTS for modelling spots are, where x is the output 
from the random number generator (0<x<l); 

(i) longitude: randomly distributed between 0° & 360° 

(ii) latitude: — § < < f , following = arcsin (2x + 1) with 
0<x<l, to eliminate an artificial concentration of spots at the pole 

(iii) spot radius: computed using the previously described log- 
normal distribution as tabulated in TableQ 

(iv) spot brightness & spot sharpness: modelled to obtain an um- 
bral to penumbral ratio of 1:3. 



2.3 Spot distributions 

The spot distributions resulting from the size distributions given in 
TableQare shown on the left hand side of FigureQ(sets 1 to 3) and 
in Figure [2] (sets 4 to 7). The degree of spottedness increases from 
the set at solar minimum (top of Figure to the set with 100% 
spot coverage (bottom of Figure |2j to show the extreme effect of 
spot coverage. In this case, as the spot coverage fraction includes 
the spot's umbra + penumbra, the result is not a completely black 
star. 



These modelled spot distributions were used to generate a syn- 
thetic lightcurve, comprising 500 points with uniform sampling and 
with random Gaussian noise (0.0004) added to match the photo- 
metric precision of our HST data. The central lightcurve wave- 
length, 4670A, was also chosen to match that of our HST data. 



3 SURFACE BRIGHTNESS IMAGE RECONSTRUCTION 

The resulting lightcurve from the modelled distributions of spots, 
as described in the previous section, is used as input to DoTS to 
reconstruct the surface brightness distribution using the maximum 
entropy x 2 minimisation method jCollier Cameronll997l) . 



3.1 System parameters 

In the reconstruction of a surface brightness distribution of a binary 
star, i t is im perative that the correct binary system parameters are 
used JVincentetalll993l) . Incorrect binary system parameters lead 
to a poor fit of the photometric lightcurve and to the appearance 
of spurious spot features on the surface brightness image to com- 
pensate for this. In this work, the assumed spot temperature used is 
1500 K cooler than the star's photospheric temperature. As the spot 
coverage increases, the surface brightness of the star will decrease 
compared to the star's photospheric brightness. This is particularly 
relevant for sets 4 to 7, where there is a high degree of spot cover- 
age. In reconstructing the surface brightness distributions for each 
of these sets the photospheric temperature, the primary radius and 
the secondary radius are treated as an unknown parameters. The 
best fitting parameters are determined using a 'grid search' method 
that uses, as input to DoTS, a grid of: photospheric temperatures 
ranging from 5500 K to 6100 K at 25 K intervals; primary radius 
values ranging from 1.22 to 1.26 Rq and secondary radius values 
ranging from 0.77 to 0.81 Rq, both in 0.005 Rq intervals. For 
each set, the lowest \ 2 value is determined by the minimum of a 
quadratic fit to these values. The resulting best-fitting parameters 
are summarised in Table [2] These values are used as the photo- 
spheric temperature and radii for the reconstruction of the primary's 
surface brightness distribution. 



3.2 Final image reconstruction 

The images reconstructed for sets 4 to 7 are shown in Figure|2| how- 
ever, the over-all spot coverage for sets 1 to 3 in FigureQwas too 
small to be reconstructed. The photometric lightcurves for sets 1 to 
3 are shown in the right-hand side of FigureQ For sets 4 to 7, the fit 
to the photometric lightcurve and reconstructed surface brightness 
images are shown respectively in the centre and the right-hand side 
of Figure|2| 
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Figure 1. (left) The distribution of spot size distributions for sets 1 to 3. (right) The photometric lightcurve resulting from these spot distributions and the Max 
Ent model fit to these lightcurves. The Observed minus Computed lightcurves are plotted at the bottom of each lightcurve. It was not possible to reconstruct 
surface brightness distributions for these three sets as the spot coverage is too low to be recovered using eclipse-mapping. 



Set 


Primary Temperature (K) 


R pr i/i?Q 


RsecAR© 


initial 


6038 


1.25 


0.81 


4 


6026± 40 


1.23()±0.003 


0.792±0.003 


5 


5948±71 


1.248±0.002 


0.804±0.003 


6 


5849± 53 


1.251 ±0.002 


0.782±0.004 


7 


5760± 45 


1.252±0.003 


0.773±0.005 



Table 2. Best fitting binary system parameters for the size distribution sets 
4 to 7 in Table 1 



What is distinctive about these plots are the spot features at 
the quadrature points (90° and 270°). These spot features are not 
artifacts of the maximum entropy reconstruction, but occur as a 



consequence of the \ 2 minimisation technique. If a star is peppered 
with small unresolvable spots, the effect on the star's lightcurve 
will be to reduce the depth of the primary eclipse. To fit this with 
\ 2 minimisation it is necessary to increase the level of the whole 
computed lightcurve to fit the reduced primary eclipse depth. The 
difference in light at the quadrature points can only be accounted 
for by placing large spurious spots at these longitudes. 



For sets 4 and 5, the spot feature at 270° longitude is fainter 
than the spot feature at 90° longitude as there is a higher degree 
of spot coverage at the first quadrature than at the second. In the 
reconstructed surface brightness distributions, the spots at 0° lon- 
gitude result from the total spot coverage between second and third 
contact, i.e. during primary eclipse, while the spots at 180° longi- 
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Figure 2. (left) The modelled distribution of spots for sets 4 (top) to 7 (bottom), (centre) The photometric lightcurve resulting from these spot distributions 
and the Max Ent model fit to these lightcurves. The Observed minus Computed lightcurves are plotted at the bottom of each lightcurve. (right) The surface 
brightness distribution that are reconstructed using the Max Ent eclipse-mapping technique, where phase runs in reverse to longitude. Note the large spurious 
spots at the quadrature points. 
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tude result from spots on the primary star visible during secondary 
eclipse. 



4 DISCUSSION 

Size distributions of starspots extrapolated from the case of the Sun 
lSolanki| [T999). have been modelled on the eclipsing binary SV 
Cam to synthesise images of stellar photospheres with high spot 
filling factors. The lightcurves resulting from the extrapolated spot 
distributions were used as input to the Max Ent eclipse-mapping 
code DoTS. The resulting surface brightness distributions, with dis- 
tinctive spot features at the quadrature points (90° and 270°), bear 
little resemblance to the input distribution of spots. These spurious 
spots at the quadrature points (90° and 270°) are not artifacts of the 
eclipse mapping technique, but a consequence of the \ 2 minimisa- 
tion method. 

The presence of spots over a large fraction of the stellar sur- 
face influences the shape of the lightcurve and consequently the de- 
rived stellar parameters. The peppering of the stellar surface with 
small unresolvable spots decreases the depth of the primary eclipse, 
resulting in a lower photospheric temperature of the primary star 
and a corresponding change in the radii of both binary components. 
The change in the binary system parameters (cf Table 2) shows how 
the degree of spottedness affects the determination of the correct 
photospheric temperature and radii when only photometric obser- 
vations are used. 

Spot features occurring at the quadrature points are commonly 
reconstructed features of RS CVn binaries and typically referred 
to as 'active longitudes'. They are interpreted as being the pre- 
ferred longitude of emergence of magnetic flux. These features 
have been observed on numerous short-period active RS CVn bi- 
nary stars suc h as EI E ri, II Pe g, a Gem, HR7275, AR Lac, SZ 
Psc, HK Lac jBerdvusina & Tuominelll998t [Lanza et aljll998t 
Lanza et all200lUjetsill99d : lHenrv et alll995tlQlah et alll99lh . 
In these cases, however, the 'active longitudes' have been deter- 
mined through a Fourier analysis of the lightcurves and clearly 
show migration over several epochs. The fixed longitudes found 
in this work would indicate that is unlikely that migration could be 
reproduced by random distributions of spots. 

Another approach for the reconstruction of a surface 
brightness distribution is the 'two-spot' model. The 'two-spot' 
model can account for the light variations of chromospherically 
active stars and exa mples of this mode l applied to SV Cam 
can b e f ound in : IZbori l & Diurasevid l2003jjA^bavral^yilJ 
1200 ll). iKiurkchieva et alJ fcOOd) iDiurasevid ll998>) " 

[ Patkos & Hempelmannl ll994l) . iBudding & Zeilikl 1 19871) and 
Zelik et al. 1 1988). However Doppler images of chromospherically 
active RS CVn stars (e.g. IPetit et alJ (2002)) show more compli- 
cated distributions of spots than can be reconstructed with the 
'two-s pot' model. To inv estigate the reliability of the 'two-spot' 
model Eat on et all 1 19961) modelled up to 40 randomly distributed 
spots on a single star, and found that they could reproduce 'two- 
spot' solutions similar to those r econstructed for actual stars. We 
have extended the work of lEaton et all 1 19961) to show that the flux 
deficit caused by the presence of many sub resolution spots can be 
modelled by two spurious spots at the quadrature points. 

It is concluded that 'two-spot' modelling or chi-squared min- 
imisation techniques are more susceptible to spurious structures be- 
ing generated by systematic errors, arising from incorrect assump- 
tions about photospheric surface brightness, than simple Fourier 
analysis of the light-curves. 
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